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Abstract 

Based on the assumption of two-quark structure of the scalar 1£q(1430), the CP-averaged 
branching ratios(BRs) and CP-violating asymmetries of charmless hadronic B q (q = u,d,s,c) —> 
ifj (1430)^0 (1430) decays are studied in the standard model(SM) by employing the per- 
turbative QCD(pQCD) factorization approach. Our predictions are the following: (1) the 
CP-averaged BRs for B q — > Kq (1430)Kq(1430) decays in both scenarios vary in the range 
of 10~ 6 ~ 10~ 4 in the SM; (2) the magnitudes of A<$ X P (B U -> i^(1430) + i?o(1430)°) and 
Aqp(B^ s —t- Kq (1430) + Xq (1430) ) in Scenario 1 are much larger than those in Scenario 2 
correspondingly; (3) there are no direct CP violations in B d)S -»■ Z^( 1430)° ^(1430) and 
B c — > Kq(1430) i^ (1430) decays in the SM because of the pure penguin and tree topology, 
respectively. A measurement of our pQCD predictions at the predicted level will favor the 
qq structure of the scalar Kq(1430) and help us to understand its physical properties and the 
involved QCD dynamics. 
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I. INTRODUCTION 



The exclusive non-leptonic weak decays of B q mesons (q = (u,d,s,c)) provide not 
only good opportunities for testing the Standard Model (SM) but also excellent places for 
probing different new physics scenarios beyond the SM. Very recently, the Belle collabo- 
ration has reported a preliminary upper limit on branching ratio of charmless hadronic 
B° /^*(1430)°Zo(1430) 01 decay (l^: 

Br(B° -> K*°K* °) = 3.21+1^+2-31 x 10 -e or < 8A x 1Q -6 at 9Q% c L ^ 

where the result is fitted for decay mode with final states K + n~ K~n + and Bt(Kq° — > 
K + rr~) w 66.7%. This measurement will be improved soon with the ongoing Large 
Hadron Collider(LHC) experiments at CERN. At LHC experiments, the b hadrons such 
as B u , B<i, B s , B c , even A& can be accessed easily. Particularly, the B c meson could be 
produced abundantly, which will make a new realm to test the SM, study the heavy flavor 
dynamics, and explore the involved perturbative and nonperturbative QCD dynamics [4j. 

In the naive quark model, Kq is a p-wave scalar(l 3 Po) particle with quantum number 
J PC = ++ . Lattice calculations (5[ on the masses of Kq and ao(1450) indicate a good 
SU(3) symmetry for the scalar sector, while the latter has been confirmed to be a qq meson 
in lattice calculations Recently, Cheng, Chua, and Yang proposed two possible 

scenarios based on the assumption of two-quark structure to describe this light scalar 
Kq in the QCD sum rule method (Iol |: the first excited state in scenario 1(S1) or the 
lowest lying state in scenario 2(S2), and made extensive studies and interesting analyses 
phenomenologically on charmless hadronic B — > Kq(P,V) (Here, P and V stand for 
the light pseudoscalar and vector mesons, respectively) decays to implicate its physical 



properties in the QCD factorization(QCDF) approach [ill, [l2j. Moreover, the people also 



made relevant investigations on Kq with other approaches/methods in hadronic B meson 



decays [13H17 



i Ply . 

181, 



It is well known that the key point of the theoretical calculation for the charmless 
hadronic B q meson decays is how to calculate the hadronic matrix element (HME) reliabl 
So far, many theoretical approaches/ methods, such as naive factorization assumption 
generalized factorization approach Jl9j], QCDF, soft-collinear effective theory(SCET) 20| , 
and perturbative QCD(pQCD) approach (21-23], are developed to make effective evalua- 
tions of HME and interpret the existing rich data. Up to now, the pQCD approach has 
become one of the most popular methods due to its unique features |24j . The annihilation 
diagrams, for example, can be evaluated here. While the strong phase for generating CP 
violation (25| in the pQCD factorization approach, is rather different from that as claimed 
in SCET H. 

Motivated by the above observations on both experiment and theory aspects, we will 
investigate the charmless hadronic B q — >■ KqK 2 decays with q — u, d, s, and c by 
employing the low energy effective Hamiltonian [23] and the perturbative QCD(pQCD) 



1 In the following section, we will adopt Kq to denote Kq(1430) + and ^0(1430)°, and K to stand for 

_ * o 

i^o(1430) and X o (1430) for convenience, respectively, unless otherwise stated. 

2 Hereafter, for the sake of simplicity, we will adopt B to denote the B u and B d mesons, unless otherwise 
stated. 
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factorization approach in this work. We here not only calculate the usual factorizable con- 
tributions, but also evaluate the nonfactorizable and the annihilation type contributions 
theoretically. We will predict the physical observables such as CP-averaged branching 
ratios(BRs) and CP-violating asymmetries in the considered decays. The large BRs and 
CP violations in the relevant considered decay channels will play an important role in ex- 
ploring the physical properties of scalar Kq. Furthermore, the pure annihilation processes 

Bd — > Kq + Kq~ and B c — > Kq + K* could provide interesting information to explore the 
underlying decay mechanism of the weak annihilation decays. 

The paper is organized as follows. In Sec. [TTJ we present the theoretical framework 
on the low energy effective Hamiltonian and formalism of the pQCD approach. Then 
we perform the perturbative calculations for the considered B q — > KqK decay channels 
with pQCD approach in Sec. IIHI The analytic formulas of the decay amplitudes for all 
the considered modes are also collected in this section. The numerical results and phe- 
nomenological analysis are given in Sec. IIVI Finally, Sec. |V] contains the main conclusions 
and a short summary. 

II. THEORETICAL FRAMEWORK 



For the considered decays, the related weak effective Hamiltonian H e g [27| can be 
written as 



#cff - ^ VQb V QD 
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+ H.c. , (2) 



*=3 

with the Fermi constant Gf = 1.16639xl0 -5 GeV -2 , Cabibbo-Kobayashi-Maskawa(CKM) 
matrix elements V, light down type quarks D = d,s, and Wilson coefficients Cj(/x) at the 
renormalization scale \i. The local four-quark operators Oi(i = 1, ■ ■ ■ , 10) are written as 

(1) Current-current (tree) operators 

O® = (D a Qp) v _A(Qpb a )v-A , = (D a Qa)v-A(Qi3bp)v-A ; (3) 

(2) QCD penguin operators 

3 = (D a b a )v~A^2(q , pq' /3 )v-A , 4 = (D a b /3 )v-A^2(^q' a ] 



V-A , 



q> q' 



(4) 



V+A , 



q' 



(3) Electroweak penguin operators 

3 3 
O7 = -{D a b a )v-A^2eq,(q'pq'p) v+A , 08 = -{D a bp) V - A e A<?Mv+A , 

3 _ q ' 3 _ q _ (5) 

°9 = ^{D a b a ) V -A^2e ql (q p q p ) v ^ A , O w = ~{D a bp) v _ A e g /(^^) y _ A . 
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with the color indices a, (3 and the notations (q'q')v±A — q'lfi(^ =t 75)9' ■ The index g' 
in the summation of the above operators runs through u, d, s, c, and b. The standard 
combinations Oj of Wilson coefficients are defined as follows, 

ai = C 2 + ^, a 2 = d + ^, ai = Ci + ^(i = 3-10) . (6) 

where the upper (lower) sign applies, when % is odd(even). Since we work in the leading 
order[C(a s )] of the pQCD approach, it is consistent to use the leading order Wilson coeffi- 



cients, although the next-to-leading order calculations already exist in the literature [27 . 
This is the consistent way to cancel the explicit /1 dependence in the theoretical formulae. 
For the renormalization group evolution of the Wilson coefficients from higher scale to 
lower scale, we use the formulas as given in Ref. [22] directly. 

The basic idea of the pQCD approach is that it takes into account the transverse 
momentum ky of the valence quarks in the calculation of the hadronic matrix elements. 
The B q meson transition form factors, and the spectator and annihilation contributions 
are then all calculable in the framework of the k^ factorization. In the pQCD approach, 
a B q — y M 2 M 3 decay amplitude is factorized into the convolution of the six-quark hard 
kernel(if), the jet function(J) and the Sudakov factor(S) with the bound-state wave 
functions($) as follows, 

A{B q M 2 M 3 ) = $ Bq ® H <g> J ® S <g) $ Ma ® ®m 3 , (7) 

The jet function J comes from the threshold resummation, which exhibits suppression in 
the small x(quark momentum fraction) region|28|. The Sudakov factor S comes from the 



k^ resummation, which exhibits suppression in the small k^ region[29|, [30(. Therefore, 
these resummation effects guarantee the removal of the endpoint singularities. 

In the practical applications to heavy B q meson decays, the decay amplitude of Eq. (jTJ) 
in the pQCD approach can be conceptually written as 3 , 

A(B q M 2 M 3 ) ~ Jd%d%d%Tr [Cit^B^h^MMQM^H^k^k^t)] ,(8) 

where fcj's are momenta of light quarks included in each mesons, and Tr denotes the trace 
over Dirac and color indices. C(t) is the Wilson coefficient which results from the radiative 
corrections at short distance. In the above convolution, C(t) includes the harder dynamics 
at larger scale than m,B q scale and describes the evolution of local 4- Fermi operators from 
uiyy (the W boson mass) down to t ~ C^^Aqcd^bJ scale, where A QCD is the hadronic 
scale. The function H(k±, k 2 , k%, t) describes the four quark operator and the spectator 
quark connected by a hard gluon whose q 2 is in the order of Aqcd^s,) an d includes 
the (^(^Aqcd^b,) hard dynamics. Therefore, this hard part H can be perturbatively 
calculated. The function $^ is the wave function which describes hadronization of the 
quark and anti-quark to the meson M, which is independent of the specific processes and 
usually determined by employing nonperturbative QCD techniques or other well measured 
processes. 



3 J(S), organizing double logarithms in the hard kernel (meson wave functions), is hidden in H (the 
three meson states). 
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Since the b quark is rather heavy, we work in the frame with the B q meson at rest 
for simplicity. Throughout this paper, we will use light-cone coordinate (P + , P~ , Pt) to 
describe the meson's momenta with the definitions P ± = Po ^ s and = (pi,p 2 )- For 

the charmless hadronic B u —y K£ K decay, for example, assuming that the Kq (K ) 
meson moves in the plus (minus) z direction carrying the momentum P2 (-P3). Then the 
two final state meson momenta can be written as 

Pi = ^§(1,1,0*.), P 2 = ^|(l-r 2 ,r 2 ,0 T ), P 3 = ^ { rll-rlO T ) } (9) 

respectively, where r 2 = m^./ms and r 3 = m-^*Jra B . Putting the (light-) quark momenta 
111 B, K* and K mesons as fei, k 2 , and k 3} respectively, we can choose 

h = { Xl P+,o,k 1T ), k 2 = (x 2 P+,o,k 2T ), k 3 = (o,x 3 P 3 -,k 3T ); (10) 

Then, for B u — > Kq + K* ° decay, the integration over k±, k 2 , and k 3 will conceptually 
lead to the decay amplitude in the pQCD approach, 

, ^0 f 

A(B U — >■ Kq K Q ) ~ / dxid X2 dx3bidbib 2 db2b 3 db 3 

xTr [C^b^x.M^k^xM 
x^ n (x 3 ,b 3 )H(x h b h t)S t {x i )e- s ^ . (11) 

where bi is the conjugate space coordinate of k^, and t is the largest energy scale in func- 
tion H(xi,bi,t). The large logarithms \n(mw/t) are included in the Wilson coefficients 



C(t). The large double logarithms (In 2 xi) are summed by the threshold resummation [28 
and they lead to the jet function S t (xi) which smears the end-point singularities on x«. The 
last term, e~ s ^\ is the Sudakov factor which suppresses the soft dynamics effectively (3lj . 
Thus it makes the perturbative calculation of the hard part H applicable at intermediate 
scale, i.e., m# scale. We will calculate analytically the function H(x{,bi,t) for the consid- 
ered decays at leading order in a s expansion and give the convoluted amplitudes in next 
section. 

In the resummation procedures, the heavy B q meson is treated as a heavy-light system( 
In the present work, the B c meson can also be viewed as a heavy-light system although c 
is the known heavy flavor quark. ). In principle there are two Lorentz structures in the B q 
meson's wave function. One should consider both of them in calculations. However, since 



the contribution induced by one Lorentz structure is numerically small [23|, [32|, [33|, and 
can be neglected approximately, we only consider the contribution from the first Lorentz 
structure 



®B q (k) 



V2Nc 



m Bq )l^(t)B q {k) 



(12) 

a/3 



where P(m) is the momentum(mass) of the B q meson, k is the momentum carried by the 
light quark in B q meson, and <f>B q is the corresponding distribution amplitude, respectively. 
In the next section, we will see that the hard part is always independent of one of the 
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K * 



V 7 

■if: m 





:f><£ I ^ 

(9) (h) 



FIG. 1: Typical Feynman diagrams for charmless hadronic B q — > KqKq decays at leading order, 
where q = u,d,s, and c, respectively. 



k + and/or k , if we make the approximations shown in the next section. The B q meson 
distribution amplitude <pB q {k) is then the function of variables A; - (or k + ) and only, 



(Ar,k 3 



d A k 

— 0b 9 (/c + , fc",k T ) 



(13) 



The transverse momentum k^ is usually conveniently converted to the b space parameter 
by Fourier transformation. 

The light-cone wave function of the light scalar Kq has been investigated in the QCD 
sum rule method as 10 



V2N C 



1Q 



* (x) + m K * (f) S K * (x) + m K * (H) - 1) 4> T K , (x) 



(14) 



a/3 



S T 

where 4>k* and 0^, , and m^* are the leading twist and twist-3 distribution amplitudes, 
and mass of the scalar Kq meson, respectively, while x denotes the momentum fraction 
carried by quark in the meson, and n = (1, 0, Or) and v = (0, 1, Or) are the dimensionless 
light-like unit vectors. 



III. PERTURB ATIVE CALCULATIONS IN THE PQCD APPROACH 

In the following, we will present the analytic factorization formulas for charmless 
hadronic B q — > KqK decays in the pQCD approach. Apart from the factorizable and 
nonfactorizable spectator diagrams, we can also calculate analytically the annihilation- 
type ones with no endpoint singularity by employing the pQCD approach. We will adopt 
F and M to stand for the contributions of factorizable and nonfactorizable diagrams from 
(V — A)(V — A) operators, F P1 and M P1 to stand for the contribution from (V — A)(V+A) 
operators, and F P2 and M P2 to stand for the contribution from (S — P)(S + P) operators 
which result from the Fierz transformation of the (V — A)(V + A) operators. 
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A. B/B s -> K* kI decays 

As illustrated in Fig. [TJ when q is it, d or s, all eight types of diagrams may contribute 
to the -B( s ) —> KqK decays. We firstly calculate the usual factorizable spectator(/s) 
diagrams (a) and (b), in which one can factor out the form factors B — » K§ and B s — > K* . 
The corresponding Feynman decay amplitudes are given as follows, 

(i) (V -A)(V- A) operators: 

/l POO 
dx x dx 3 J M&1M&3 

x0 B(s) (xi, 61) { [(1 + x 3 )0s(x 3 ) + r 5 (l - 2x 3 )(0|(x 3 ) + <f>s(x 3 ))] 
xh fs (x 1 ,x 3 ,bi,b 3 )E fs (t a ) + 2r s 4>i(x 3 ) h fs (x 3 ,x 1 ,b 3 ,b 1 ) E fs (t b )} ,(15) 

(ii) (V -A)(V + A) operators: 

Ff s l = F fs , (16) 

(iii) (S — P){S + P) operators: 

r 1 poo 

Ff s 2 = 16irC F f s m 2 B r s / dx x dx 3 I b x dbib 3 db 3 

Jo Jo 

x<p B(s] (xi, 61) { [(f> s (x 3 ) + r s [(2 + x 3 )(f>s(x 3 ) - x 3 f(x 3 )]] 
xh fs (x 1 ,x 3 ,b 1 ,b 3 )E fs (t a ) + 2r s <j%(x 3 ) h fs (x 3 , x 1 ,b 3 ,b 1 ) E fs (t b )} ;(17) 

where rs = ms/ms (Hereafter, for simplicity, we will use S to denote K$ and its charge 
conjugation K* in the explicit expressions of factorization formulas.) and Cf = 4/3 is a 
color factor. The convolution functions Ei, the factorization hard scales ti, and the hard 



functions hi can be referred to Ref. [34j. 

For the non-factorizable spectator (nfs) diagrams 1(c) and 1(d), the corresponding 
decay amplitudes can be written as 

(i) (V -A)(V- A) operators: 

32 f 1 f°° 

M n f s = —tcC fvt?q \ dx\dx 2 dx 3 \ b\db\b 2 db 2 

V6 (s) Jo Jo 

x0b (s) (a?i, b 1 )(t) S (x 2 ) { [(1 - x 2 )(j)s(x 3 ) - r s x 3 ((j) s s (x 3 ) - ($(x 3 ))] 

xE nfs (t c )h c nfs (x 1 ,x 2 , x 3 , bi, b 2 ) - hf tfs (x 1 ,x 2 , x 3 , bi, b 2 ) 

x [{x 2 + x 3 )0s(x 3 ) - r s x 3 {<j) s s (x 3 ) + f s (x 3 ))] E nfs {t d )} , (18) 
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(ii) (V - A)(V + A) operators: 

32 f 1 f°° 

M^} s = —=ixCFm 2 B j dx\dx 2 dx 3 I bidbib 2 db 2 
V6 <s) Jo ' Jo 

x0 B(s) (xi, h)r s { [(1 - x 2 )(i)s(x3)((i)s(x2) + <PsM) 

-r s (<j4(x 2 )[(x 2 -x 3 - l)<l4M ~ (x 2 + x 3 - l)(j) T s {x 3 )\ 

+f s {x 2 )[{x 2 + x 3 - l)0|(x 3 ) + (1 - x 2 + x 3 )<J>s])] E nfs {t c ) 

xh° nfs (x 1 ,x 2 , x 3 , bx, b 2 ) - h* fs (x 1 ,x 2 , x 3 , b x , b 2 )E nfs (t d ) 

x [x2((p S si,x 2 ) - (f>l(x 2 ))(j)s(x 3 ) + r s (x 2 ((f) s s (x 2 ) - <t> T s (x 2 )) 

x (0f (x 3 ) - 0l(x 3 )) + x 3 (0| (x 2 ) + 0l(x 2 ))(0|(x 3 ) + 0l(x 3 )))] } ,(19) 

(iii) (S — P)(S + P) operators: 

32 Z" 00 
M^ 2 S = —nC F m 2 B / dx\dx 2 dx 3 I hdbib 2 db 2 (f)B (s) (xi,h) 

x<Ps(x 2 ) { [(x 2 -x 3 - i)<f> s (x 3 ) + r s x 3 ((p s s (x 3 ) + <£s(z 3 ))] 

xE nfs (t c )hl fs (xi,x 2 , x 3 , b 1: b 2 ) + h d nfs (xi, x 2 , x 3 , h, b 2 ) 

x [WsrfoO - r s x 3 (^f (x 3 ) - (f>s(x 3 ))] E nfs (t d )} ; (20) 

In the above three formulas, i.e., Eqs. (fT8l) - (l20l) . one can find that there exist cancelations 
between the contributions of the two diagrams in Fig. [T](c) and[T](d). 

The Feynman diagrams shown in Fig. [T](e) and []Jf) are the non-factorizable 
annihilation(n/a) ones, whose contributions are 

(i) (V -A)(V- A) operators: 

32 f 1 f°° 

M nfa = -=rcC F m 2 B / dx x dx 2 dx 3 \ b 1 db 1 b 2 db 2 

V6 {s) Jo Jo 

x< Pb {s) (xi, b x ) { [(1 - x 3 )(j) S (x 2 )(j) S (x 3 ) - r s r s (</>|(x 2 ) 
x [(1 + x 2 - x 3 )4>s(x 3 ) - (1 - x 2 - x 3 )(pl(x 3 )] + 4> T s {x 2 ) 
x [(1 -x 2 - x 3 )(j) s s {x 3 ) - (1 + x 2 - xa)4>s(x 3 )])] E nfa (t e ) 
xhl fa (x u x 2 , x 3 , bl , b 2 ) - E nfa {t f )h{ fa { Xl ,x 2 , x 3 , bl , b 2 ) 
x [x 2 (j) S (x 2 )(j)s(x 3 ) + r s r s (jc 2 )[(a? 3 - x 2 - 3)0§ (x 3 ) 
-(1 -x 2 - x 3 )(j} T s {x 3 )} + (j) T s {x 2 )[(l -x 2 - x 3 )(j) s s {x 3 ) 
-(l-x 2 + x 8 )^(x s )])]} , (21) 
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(ii) (V -A)(V + A) operators: 

32 f 1 f°° 

M n f a = —^irCpin^ I dx\dx2<lxz l bidbib2<ib2 
v6 s) Jo Jo 

X(f) B(s) (x 1 ,b 1 ) { [r s x 2 (i)s(x 3 )((f) s s (x2) + <Ps(x 2 )) + r 5 (l - x 3 ) 
x0 5 (x 2 )(0|(x 3 ) -(f>s(x 3 ))] E nfa (t e )h e nfa (x 1 ,X2 } x 3 ,b ll b 2 ) 
+ [r s (2 - x 2 )((f> s s (x 2 ) + 0£(x 2 ))0 5 (x 3 ) + r 5 (l + x 3 ) 

xfaMttsM -4>s( x 3))] E nfa (t f )h f nfa (x 1 ,X2,x 3 ,b 1 ,b 2 )} , (22) 

(iii) (S — P)(S + P) operators: 

32 f 1 f°° 

Mnfa = —/=^CFm 2 B / dx±dx 2 dx 3 / bidbib 2 db 2 
v6 <a) Jo Jo 

x0£t (s) (xiA) {[(1 - ^3)0sO2)0sO 3 ) -r 5 rs (0f(z 2 ) 

x [(x 2 - x 3 + 3)0f (x 3 ) - (1 - ar 2 - 2; 3 )0s(x 3 )] + 0|(x 2 ) 

x [(1 - x 2 - x 3 )4> s s (x 3 ) + (1 - x 2 + x 3 )</>g(x 3 )])] E nfa (t f ) 

Xhi fa ( Xl , X 2 , X 3 , h, b 2 ) - E nfa (te)K fa (x 1 ,X2, X 3 , b U b 2 ) 

x [x2(p s (x 2 )(ps(x 3 ) + r s r s (</>f {x 2 )[{x 3 - x 2 - l)^f (a: 3 ) 

-(1 - £ 2 - a; 3 )05(x 3 )] + 0f (x 2 )[(l -x 2 - x 3 )(j)§ (x 3 ) 

+ (l + x 2 -x 3 )0i(x 3 )])]} ; (23) 

For the last two diagrams in Fig. (TJ i.e., the factorizable annihilation(/a) diagrams 
UX.&) and [0(h), we have 

(i) (V - A)(V - A) operators: 

/l poo 
dx 2 dx 3 J b 2 db 2 b 3 db 3 {[x 2 5 (x 2 )0 5 (x 3 ) - 2r s r s 

x ({x 2 + l)0f (ar 2 ) + (ar 2 - l)0|(x 2 )) #f (ar 3 )] /i /a (x 2 , 1 - z 3 , 6 2 , 6 3 ) 
xE fa {t g ) - [(1 - x 3 )0 5 (x 2 )0 s (x 3 ) + 2r s r 5 0f (x 2 ) ((x 3 - 2)0f (x 3 ) 
-x 3 0g(x 3 ))] Ef a (t h )hfa{l - x 3 ,x 2 ,b 3 ,b 2 )} , (24) 

(ii) (V - + A) operators: 

Ffa 1 = F fa , (25) 

(iii) (S — P){S + P) operators: 

F^ 2 = -167rC F m| / dx 2 dx 3 / b 2 db 2 b 3 db 3 { [2r s 5 (x 2 )0f (x 3 ) 

Jo Jo 

-r-sfX 2 (0f (x 2 ) - (pl(x 2 ))(l)s(x 3 )] h fa (x 2 , 1 - x 3 , b 2 , b 3 )E fa (t g ) 
+ [r s (l - x 3 )0 s (x 2 )(0|(x 3 ) + 0i(x 3 )) - 2r s 0f (ar a )^(ar 3 )] 
xE fa (t h )h fa (l - x 3 ,x 2 ,b 3 ,b 2 )} . (26) 



It is interesting to notice that there is a large cancelation in the factorizable annihilation 
Ff a , i.e., Eq. (|24p . from the diagrams QJg) and[TJh), which can result in the tiny or small 
deviations from zero and can be seen numerically as displayed in the last column of the 
2nd line of Tables HI and ITT1 In the SU(3) limit, in particular, this cancelation will lead to 
the exact zero contribution. 



B. B c K^ + K* ° decay 

In the SM, charmless hadronic B c — > Kq K decay can only occur through the 
pure annihilation- type diagrams. From the effective Hamiltonian fl2]), there are 4 types 

,,,0 

of diagrams contributing to the B c — > Kq K decay as illustrated in Fig. []7e)-(h), 

which result in the Feynman decay amplitudes Fj a and M' n j a from singly current-current 

operators, respectively. Following the same procedure as stated in the above subsection, 

*o 

we can obtain the analytic decay amplitudes for B c — > Kq K mode, 

/»1 poo 

F'f a = — 8irC fttVq^ / dx 2 dx 3 / b 2 db 2 b 3 db 3 
c Jo Jo 

x {hf a (l - x 3 , x 2 , b 3 , b 2 )E' }a {t' g ) [x 2 <Ps O2) <Ps(x 3 ) - 2r 5 r 5 0|(x 3 ) 

x ({x 2 + 1)<P S S {x 2 ) + (x 2 - l)0l(x 2 ))] - h fa {x 2 , 1 - x 3 , b 2 , b 3 )E' fa (t' h ) 

x [(1 - x 3 )<f> s (x 2 )<t> s {x3) + 2r 5 r 5 0|(x 2 ) ((x 3 - 2)</»|(x 3 ) - x 3 ^{x 3 ))] } ,(27) 

Mnf a = TrCFfTi 2 Bc / dx 2 dx 3 / b\dbib 2 db 2 

x {h'lfai^ x 3 , b u b 2 )E' nfa (t' e ) [{r c - x 3 + l)4>s{x 2 )(j)s{x 3 ) - r s r s (</>f(z 2 ) 
x ((3r c + x 2 - x 3 + l)0f (x 3 ) - (r c - x 2 - x 3 + 1)4> T S (x 3 )) + 4> T s {x 2 ) 
x((r c -x 2 -x 3 + l)0f(x 3 ) + (r c -x 2 + x 3 - l)<f%(x 3 )))] - Kfatff) 
x [{r b + r c + x 2 - l)0 s (x 2 )0 5 (x 3 ) - r s r s (0f (x 2 )((4r b + r c + x 2 - x 3 
-l)(f>s(x 3 ) - (r c + x 2 + x 3 - l)0|(x 3 )) + 05(x 2 )((r c + x 2 + x 3 - 1) 
x0|(x 3 ) - (r c + x 2 -x 3 - l)0f(x 3 )))] /i'i /a (x 2 ,x 3 ,6i,6 2 )) , (28) 



where the non-relativistic approximation form of the distribution amplitude <pB c for B c 
meson has been used, the convolution factor E[, the hard scale t' i: and the hard function 



hj are referred to Refs. [17|, |35j. Moreover, r b = rrih/mB,., r c = m c /rriB c , and r b + r c « 1 
in 5 r meson. 



C. Decay Amplitudes for i? g — >■ KQK Q (q = u, d, s, c) Channels 

By combining various of contributions from the relevant Feynman diagrams together, 
the total decay amplitudes for the four penguin-dominated decays — > KqK and the 

pure annihilation processes B d — > Kq + Kq~ and B c — > Kq + K* ° can then read as, 
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1. The total decay amplitudes of B — > KqK* decays: 



A(B U -+ K* + K ) = \ u 



M nfa d 



- A f 



F fs (a 4 - iaio) + F/^K - \a%) 



+M nfs (C 3 - ic 9 ) + M nfs (C 5 - ic 7 ) + M nfa 
x (C 3 + C 9 ) + M# a (C 5 + C 7 ) + f B Ff a 2 {a 6 + a 8 ) 



, (29) 



where A„ = V* b V ud and At = V* h V tf i, 



t = v t b v td- 



A{B d tf * + tf *-) = A, 



M nfa C 2 



At 



M n/a (C 4 + C 10 ) + M„ p / a (C 6 + C 8 ) 



+M n/a [^* + X*-](C 4 - -C 10 ) 



(30) 



„4(5 d ) = -At 



Ffs{a A - -aio) + Ff s 2 (a 6 - -a 8 ) + (M n/s + M n/o ) 



x (C7 S - ±C 9 ) + (M^ + M„7J(C 5 - l -C 7 ) + (M nfa 



1 



+ [K* U Kl"])(C 4 - ^C 10 ) + {Mnf a + W » Kll) 



-*o n 



x(C 6 -^ 8 ) + / B F; a 2 (a 6 -i C 



(31) 



2. The total decay amplitudes of B s — > KqK* decays: 



A{B S K * + K *-) = A; 



F /S ai + M nfs d + M nfa C 2 



-a; 



F fs (a 4 + aio) 



+Ff/(a 6 + o 8 ) + M„ /s (C 3 + C 9 ) + M^(C 5 + C 7 ) 
+M n/a (C 3 - l -C, + C A - ho w ) + M nfa [K* + K*-} 



x (C 4 + C w ) + M# (C 6 - ^-C 7 ) + M# a (C7 6 - 



1„, 1 



+M n p / a [K* + K* -]{C 6 + C 8 ) + /B s F/ a 2 («e - §a 8 ) 



P2/ 



(32) 
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where = V^ b V us and X' t = V t * b V ts . 

-K 



A(B S -> K*°Kl°) 



F fs (a 4 - iaio) + Ff s 2 (a 6 - ]-a 8 ) + (M nfs + M nfa ) 



2 «./ /s v u 2 
x(C 3 - -C 9 ) + (M„ p / S + M n 7J(C 5 - ic 7 ) + (M n/a 
+[< o <])(C 4 - ~C 10 ) + (M„ p / a + [< <]) 



x(C 6 -^ 8 ) + / Bs F; a 2 (a 6 -ia 8 ) 



(33) 



*o 



3. The total decay amplitude of B c — > Kq K decay: 



A{B C 



Ok 



fB c F'f a ai + M' nSa C x 



(34) 



In the above decay amplitudes for the channels B^a — >■ KqK q , based on the relevant 
discussions after Eq. f l2"B"j) . we have neglected the factorizable annihilation contributions 
Ff a in Eqs. ( l2T?j) -f l3"3"j) induced from the small SU(3) symmetry breaking effects. However, 

for the pure annihilation mode B c — > Kq K , to present the large annihilation contribu- 
tion occurred in this considered B c decay channel, we therefore have kept the factorizable 
decay amplitude F'? a in Eq. (134"1) . 



IV. NUMERICAL RESULTS AND DISCUSSIONS 

In this section, we will make the theoretical predictions on the CP-averaged BRs and 
CP- violating asymmetries for those considered B q — > KqK decay modes. In numerical 
calculations, central values of the input parameters will be used implicitly unless otherwise 
stated. Firstly, we shall make several essential discussions on the input quantities. 



A. Input Quantities 



The pQCD predictions depend on the inputs for the nonperturbative parameters such 
as decay constants and universal distribution amplitudes for heavy pseudoscalar B q and 
light scalar Kq mesons. For the B/B s mesons, the distribution amplitudes in the b space 
have been proposed 



4>B(x,b) = Nbx 2 (1 — x) 2 exp 



in Refs. [21], |22j and 



N Bs x 2 (l — x) 2 exp 



V 



xrriB 



UJb J 



xm Bs 

Ubs 



u 2 b 2 



(35) 



colb 2 



(36) 
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in Ref. [33J], respectively, where the normalization factors Nb (s) are related to the decay 
constants , through 



B 



dxcf> B (x,b = 0) = ^§ . (37) 



o 



2^6 ' 



In recent " lot of studies for B decays have been performed in the pQCD ap- 



proach [2l|, |22j, and the shape parameter u b has been fixed at 0.40 GeV by using the 
rich experimental data on the B mesons with /# = 0.19 GeV. Correspondingly, the nor- 
malization constant Nb is 91.745. For B s meson, considering a small SU(3) symmetry 
breaking, since s quark is heavier than the u or d quark, the momentum fraction of s 
quark should be a little larger than that of u or d quark in the B mesons, we therefore 
adopt the shape parameter co bs = 0.50 GeV [33] with f Bs = 0.23 GeV, then the corre- 
sponding normalization constant is Nb s = 63.67. In order to analyze the uncertainties of 
theoretical predictions induced by the inputs, we can vary the shape parameters u b and 
u bs by 10%, i.e., u b = 0.40 ± 0.04 GeV and u bs = 0.50 ± 0.05 GeV, respectively. 

As for the double heavy-flavored B c meson, since it embraces two heavy quarks b and c 
simultaneously, the distribution amplitude 0b c would be close to 5(x — m c /m Bc ) [36| in the 



non-relativistic limit, we therefore adopt the non-relativistic approximation form(See [35 
and references therein), 

<t>B c { x ) = ^% S ( x - m c/ m B c ) , (38) 

where fs c is the decay constant for B c meson and to be determined by the precision data 
in principle. Unfortunately, however, there are no available experimental measurements 
on fs c , we then adopt the result calculated in Lattice QCD [37], 

f Bc = (489±4±3) MeV. (39) 

For the scalar Kq, its leading twist light-cone distribution amplitude <f>K*(x,fi) can be 
generally expanded as the Gegenbauer polynomials [13, 38 ] : 



<f>K S (x,») = ^x(l-x)|/^(//)+/^(/i)£s m (//)^/ 2 (2x-l)}, (40) 

where fx*{^) and /jf*(//), B m (fi), and Cm 2 (t) are the vector and scalar decay constants, 
Gegenbauer moments, and Gegenbauer polynomials, respectively. 

There exists a relation between the vector and scalar decay constants, 

itlk* 

fx* = Vk*Jk* and [i K * = — -- — °- — -- , (41) 

where mi and m 2 are the running current quark masses in the scalar Kq. 

The values for scalar decay constants and Gegenbauer moments in the distribution 



amplitudes of Kq have been investigated at scale fi = 1 GeV in scenarios SI and S2 [10 



f K « = -0.300 ± 0.030 GeV, B x = 0.58 ± 0.07, B 3 = -1.20 ± 0.08 (SI), 
f K , = 0.445 ± 0.050 GeV, B 1 = -0.57 ± 0.13, B 3 = -0.42 ± 0.22 (S2) ; (42) 
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As for the twist-3 distribution amplitudes (p s K * and 4>k*i they have been investigated 
in only S2 with large uncertainties [39]. We therefore adopt the asymptotic forms in our 
numerical calculations for simplicity: 



■Ik*-. 



1 



U = — 7 =f K *(l - 2x). 



(43) 



The QCD scale (GeV), masses (GeV), and B q meson lifetime(ps) are [2ll. |22|. |40 
A. 



MS 

m B 



0.250 
5.279 
1.638 



m w 
m b 

T B d 



80.41 
4.8, 
1.53 , 



6.286 , m Bs = 5.366 ; 
1.5 , m K . = 1.425 ; 

1.47 , r Bc = 0.46 . (44) 



For the CKM matrix elements, we adopt the Wolfenstein parametrization and the 
updated parameters A = 0.814, A = 0.2257, p = 0.135, and fj = 0.349 . 

Utilizing the above chosen distribution amplitudes and the central values of the relevant 
input parameters, we can get the numerical results in the pQCD approach for the form 

factors Fq ^ K ° and Fq[~* 1x 



from Eq. (fT5l) at maximal recoil as follows, 



<r*v=o) 



n Q/l+0.04/ . \+0.03/f N+O.Ol/oS 1 



0.63 



1-0.10/, , V 

5 7+ 09 ( 



-0.07? j \+0.06/ R S , ^ 
0.07 U5j -0.06 



1+0-03/ 7 \ 
I-0.03USJ 
^+0.06/ J > 
)-0.07WS; 



-0.01 V^i 
f0.05/ r>S 

~omy JD i 



(SI) 
(S2) 

(SI) 
(S2) 



(45) 
(46) 



where the errors arise from the shape parameter Ub(ui, s ) in B(B S ) meson distribution 
amplitude, the scalar decay constant f$, and the Gegenbauer moments Bf in the light 
Kq distribution amplitude, respectively. These values agree well with those as given in 
Ref. & 



B. CP-averaged Branching Ratios 

In this subsection, we will analyze the CP-averaged BRs of the cahrmless hadronic 
B q KqKq decays in the pQCD approach. For B q — » KqK q decays, the decay rate can 
be written as 

r = ^^(1 - 2r 2 s )\A(B q -+ K*K* )\ 2 , (47) 

where the corresponding decay amplitudes A have been given explicitly in Eqs. (I29ti34j) . 
Using the decay amplitudes obtained in last section, it is straightforward to calculate 
the CP-averaged BRs with uncertainties as displayed in Eqs. (Fi8" ]) - (j5T?]) . AST]) - (152")) . and 
(I53p . The dominant errors are induced by the uncertainties of the shape parameters 
Ub = 0.40 ± 0.04 GeV for B mesons and Ub s = 0.50 ± 0.05 GeV for B s meson, the charm 
quark mass m c = 1.5 ± 0.15 GeV for B c meson, the scalar decay constants fs and the 
Gegenbauer moments Bf(i = 1,3) for the scalar Kq, and CKM matrix elements (p,fj), 
respectively. It is worth of mentioning that the variation of the CKM parameters has a 
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little effects on the BRs of these considered B q — > KqK* decays in the pQCD approach 
and thus will be neglected in the numerical results as shown in Eqs. (I48p - (l53p . 

The pQCD predictions for the CP-averaged BRs of the decays under consideration 
within errors in both scenarios SI and S2 are the following, 



Br(B _> k*+k*°) ~ J i-iioiH)io:4(/5)^i(sf) x 10- 5 (si) 

Br[B u -+ K A ) ~ | 1.9i°;Ko; 6 )i^(/5)±^(5f) x 10" 5 (S2) 
^ / l.llg:i(a; 6 )l°;i(/ s )^(Sf) x 10- 5 (SI) 



( 4 * 



Br(S„ Jt; K ) » <| 2 : i|?(a; 6 )ii°(/ s )+^(Sf) x 10" 5 (S2) ' (49) 
Z / U ,/ A " /V ° ; ~ 1 2.2^K)^i(/ 5 )^;«(5f) x 10- 6 (S2) ' (50) 
. „- <■ ^ / 2.5l°j(o; 6s )^i(/ 5 )i°i(Sf) x 10- 4 (SI) 



Br{B s ^K K ) „ <; 5 .4i°i(o; 6s )l|f (/ 5 )^|(5f) x 10~ 4 (S2) ' (51) 



rwr ~ / 2 - 3 -oi( w ^)^oi(/<s')^o:4(^ 5 ) x 10 4 (SI) 

£r(^ A A ) ~ I 5.2i°J(a; 6s )l 2 2; «(/ 5 )^(5f ) x 10" 4 (S2) 

*<S „ / 2.HL 7 (^)lo:?(/5)iS:I(5f) x 10- 4 (si) 



(52) 



A K ) ~ I 3.0i°; 6 7 (m c )l^(/ 5 )l?|(5f ) X KT 5 (S2) " (53) 

It is easy to see that above BRs are rather large, in the range of to 10~ 4 , and can 
be measured at B factories for the relevant B Ujd decays, and at the LHC experiments 
for the considered B s /B c decay modes. More importantly, the pure annihilation decays 

B d -> K* + K*~ and B c ->■ A* + A with large BRs can provide important information 
about the weak annihilation amplitudes and more evidences on the sizable annihilation 
contributions in B physics, then further shed light on the corresponding annihilation decay 
mechanism. 

By comparison with those of B/B s — > KK decays jH, 4l|- 44 ] , one can easily find that 



the pQCD predictions for the BRs of B/B s — > KqK channels are much larger than that 
for the corresponding B/B s —» KK decays by generally a factor of 10. For the pure 
annihilation mode B d — > Aq + Aq~, however, the enhancement factor is about 100. The 
main reason for so large difference is that the QCD behavior for the p-wave scalar Kq 
is very different from that for the s-wave pseudoscalar kaon, which can be seen clearly 
from their distribution amplitudes (lol . 45]: the former governed by the odd Gegenbauer 



moments, while the latter dominated by the even ones, apart from the small symmetry 
breaking term af. 

As mentioned in the introduction, up to now, there is only one preliminary upper limit 
on the branching ratio of B° — > Kq°K* at the 90% confidence level 



Br{B° -> K*°K* °) < 8.4 x 10" 6 , (54) 



While in the pQCD approach, our theoretical predictions of the BRs for B° — > Kq°K* ° 
decay in both scenarios within errors are as follows, 

Br(B° -> K*°K*°) ~ I 6 ~ 18 X 10_6 (S1) (55) 
BryB — > A A j ~ 1 6 ~ 40 x lO" 6 (S2)' 
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It is easy to see that the pQCD predictions within theoretical errors in both scenarios 
are consistent with currently available experiment upper limit, and will be tested directly 
when more data samples are collected at the LHC experiments. 

As shown in Eqs. (|48l) - (l53l) . the pQCD predictions have a strong dependence on the 
input parameters describing the nonperturbative behavior of the light scalar Kq ( the 
scalar decay constant fs and Gegenbauer moments Bf), and a moderate dependence on 
the shape parameters Ub(uib s ) in and the charm quark mass m c in <p Bc . Therefore, 

once the relevant experiments could provide the precise measurements for B — > KqK 
modes, one can has a better understanding about the nonperturbative hadron dynamics 
mK*. 

In order to facilitate the discussion, we present the central values of our predictions 
in terms of the topological amplitudes in Tables U and HI1 where J-/ s , J-/ a , M. n fs and 
M-nfa denote the decay amplitudes from factorizable spectator, factorizable annihilation, 
nonfactorizable spectator, and nonfactorizable annihilation contributions, respectively. 

TABLE I: The factorization decay amplitudes (in unit of 10 -3 GeV 3 ) of the charmless hadronic 
B q — > KqKq decays in SI, where only the central values are quoted for clarification. 



Channels 




M nfs 


M^\ 

nfa 




B u -> K* +Kl U 
B d -> K* + K*~ 
B d -)> KfKf 


1.246 - i0.518 
0.0 

1.246 - i0.518 


0.468 - i0.729 
0.0 

0.468 - i0.729 


-0.733 + i0.335 
2.070 - «2.517 
0.269 - il.118 


2.540 + *1.614 
-0.013 + i0.009 
2.690 + il.741 


B s -> K*+K*- 
B s -+ KfKf 


-6.020 + i0.190 
-6.270 


-0.736 + il. 594 
-2.732 + i2.310 


-2.267 + »5. 741 
-2.727 + i6.376 


-10.190 - il4.192 
-10.187- il4.273 


B c -+ K*+K* U 


0.0 


0.0 


35.744 - ill. 105 


1.658 - i2.530 



TABLE II: Same as Table H but in S2. 



Channels 




Mnfs 


M<\ 

nfa 




B u -)> K* + K*v 
B d -> K* + K* - 
B d KfK* Q ° 


2.759 - zl.147 
0.0 

2.759 - il.147 


-0.364 + i0.377 
0.0 

-0.364 + i0.377 


0.601 + i0.407 
-1.934 + i0.527 
-0.800 + il.260 


2.528 + ?4.018 
0.004 - iO.001 
2.629 + i4.035 


B s -+ K* + K* Q - 
B s -)> KfKf 


-13.606 + i0.506 
-14.158 


0.833 - il.071 
2.289 - i0.961 


4.562 - i4.863 
4.968 - z4.847 


-5.022 - i25.089 
-5.293 - i25.192 


B c -+ K* + Kl V 


0.0 


0.0 


-6.626 - il4.063 


1.387 + »0.083 



Based on the numerical results as shown in Tables [I] and [III one can easily observe that 
the factorizable contributions Tf s and Tf a govern the considered four penguin-dominated 

decay modes B u -+ K* + K* ° ', B d -+ K*°K* ° ', B s -+ K* + K*-, and B s -+ K*°K* ° ', where 



16 



both J-// and Tf a are mainly determined by the contributions induced by (S — P)(S + P) 
operators, e.g., see Eqs. (1291) and ( I3TH331) : while the nonfactorizable contributions -M^L 

dominate the rest pure annihilation channel B d — > Kq + Kq~ (B c — > Kq + K* °). The small 
magnitudes of the factorizable annihilation contributions in B d — > Kq + K^~ and B c — > 

KqKq decays are induced by the tiny SU(3) symmetry breaking effects. The difference 
between Ff a (B d — > Kq + Kq~) and T'^ a {B c — > Kq + K* ) in two scenarios as presented in 
Tables HI and HT] however, is mainly governed by two factors apart from /b c //b ~ 2.6: 
one is the ratio of V c b/V u b ~ 11.5, the other is the large ratio of a\/a2 ~ 1000 at the 
characteristic hard scale m^jl [33j j. The above phenomenological analysis can also be 
applied to the nonfactorizable annihilation topology. Anyway, one can see that these two 
considered pure annihilation modes are determined by the contributions arising from the 
nonfactorizable diagrams, i.e., M. n fa and Ai' n f a (See Tables HI and HT1) . respectively. 

It should be stressed that all these considered B q — > KqK decays have large weak 
annihilation contributions, which can be seen clearly in Tables H] and [Til It is expected 
that these decay channels with precision measurements at the experiments will offer a 
good platform to study the underlying annihilation mechanism. 

Using the pQCD predictions, we get the following interesting information on the ratio 
of the BR's for two sets of considered decays in both scenarios SI and S2, 

TB* BrjB u -> if * + ^°) pQCD „ Br(B s -> K*°K* \qcb „ 1 ( . 

Tb u Br{B d -+ KfT Q \ QCD ~ ' ' Br(B s -+ ^+^-) p qcd ~ 

If these ratios are measured in the near future, they can offer a great opportunity to study 
the QCD dynamics involved in these four decay channels and will be helpful to test the 
adopted scenario of the scalar Kq in this work. 

On the other hand, it is very interesting to notice the ratios among the theoretical BRs 
of the considered modes in the pQCD approach 



Br{B u -> K* 




)S2 


Br{B u -> K* Q 


^0 


)si 


Br(B s K* 


^0 


)S2 


Br(B s -> K* 




)si 


Br(B d -> K* + 


K*-) 


SI 


Br(B d -> K* + 




S2 



Br(B d -> KfK* ) S2 



Br(B d — >■ Kn°K, 



ii /v o°)si 



Br{B s #S°# Q 



;s2 



(1.7 -1.8), (57) 



(2.2 ~ 2.3) , (58) 



Br(£ c -> K* + irn°)si , s 

1.5 , — — _y ps 7.0 . (59) 

Br(B c -> K * + #* ) S2 

where the central values for the CP-averaged BRs have been quoted. As given in Eqs. ( 1571) - 
(1591) . one can easily find that the CP-averaged BRs for the weak annihilation processes 



4 For the considered decays, with scalar decay constant fs contributes to the magnitude of Ff s 
dominantly. This is just because the contributions arising from the vector current operators are very 
small with the strongly suppressed vector decay constant fs- 
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B d -> K* + K*- and 5 C ->■ K* + K in SI are larger than those in S2 to different extent. 

While for other four penguin- dominated B u -> K* + K , fi d K*°K* , 5 S -)• #5°#J , 
and B s — >• A^q + A"q~ channels, the CP-averaged BRs in the second scenario are larger than 
those in the first one with a factor around 2. These two different patterns might reveal 
the different QCD dynamics involved in the corresponding decay channels. The above 
relevant relations can be confronted with the near future experiments. 



C. CP- violating Asymmetries 

Now we turn to the evaluations of the CP- violating asymmetries for B q — > K^K* 
decays in the pQCD approach. 

For the charged B u and B c decays, the direct CP-violating asymmetry A^p can be 
defined as: 

Acp ~ \M\ 2 + \M? ' { } 



where M denotes the decay amplitude of charged -B u ( c ) — > Kq K decays, while M 
stands for the charge conjugation one correspondingly. 

Using Eq. fIBTJ]) . it is easy to calculate the direct CP- violating asymmetries as listed in 
Eq. (IBT]) for the considered B u decay in SI and S2, 

CP{ U ° ° j ~ I — 24.4±JJ (u b ) to'Kfs) tti (Bf) t\i (CKM) % (S2) ^ 

The large CP-violating asymmetries flBTj) plus large branching ratios ( 148)) in both scenarios 
are clearly measurable in the B factories and LHC experiments. If these physical quantities 
could be tested at the predicted level, it is doubtless that one can determine the better 
scenario of Kq meson and further understand the involved QCD dynamics. 

Because only the tree topology is involved, there is no direct CP violation in the 

considered B c decay mode, i.e., A C J P {B C — > Kq + K* ) = in both scenarios. 

As for the CP- violating asymmetries for the neutral decays Bd( s ) — > KqK , the effects of 
B d ^—B d ^ mixing should be considered. The CP-violating asymmetries of B d r s }(B d ( s )) — > 

Kq + Kq~, Kq°K* decays are time dependent and can be defined as 

A _ r (B d{s) (At) -> f CP ) - r {B d(s) (At) -> fcp) 
CP ~ r (B d{s) (At) f C p) + r (B d{s) (At) fcp) 

= Acp cos(Am (s) At) + A^p sin(Am (s) At), (62) 
where Am,( s ) is the mass difference between the two -S°(s) mass eigenstates, At = t C p — t tag 

is the time difference between the tagged -B^(s) (B d ^) and the accompanying B d ^ (B^,*.) 
with opposite b flavor decaying to the final CP-eigenstate fcp at the time tcp. The 
direct and mixing induced CP-violating asymmetries A^piCf) (or Af in term of Belle 
Collaboration) and A^^Sf) can be written as 

A dir n l^cpf ~ 1 , mix c 2Im(A CP ) 

A CP - L f - i I I \ 12' A cp - - i , I \ |i; l bd J 

1 + | Acp J- + | Acp 



18 



with the CP-violating parameter \qp 



Vtb V td(s) (fcp\H ef f\B d{s] ) 



CP = Vf v v* ' 77 — i77 — Tr^T' ^ ^ 

where T]f is the CP-eigenvalue of the final states. Moreover, for 5 S meson decays, a non- 
zero ratio (Ar/r) Ss is expected in the SM 46_|, |47|. For £> s — y KqK* decays, the third 



term *4.Ar s related to the presence of a non-negligible Ar s to describe the CP violation 
can be defined as follows |47j : 

^ = ^-1, (65 ) 



1 + A 



CP 1 



The three quantities describing the CP violation in B s meson decays shown in Eqs. (I63I) 
and (165]) satisfy the following relation, 

\A% P \ 2 + \A^\ 2 + \A Ar f = 1. (66) 

For B9sfl^ 8 ) —t- Kq°K* ° decays, they do not exhibit CP violation in both scenarios, 
since they involve only penguin contributions at the leading order in the SM, as can be 
seen from the decay amplitudes as given in Eqs. d3I]) and (J33D- Then for B°JB S K*°K* 
mode, according to Eq. ( )66|) . the third term of the CP violation Aat 3 = 100%. If the 

experimental data for the direct CP asymmetries Aq P in B d ^ — y Kq°K* decays exhibit 
obviously nonzero, which will indicate the existence of new physics beyond the SM and 
will provide a very promising place to look for this exotic effect. 

For B® s j/B( s j — y Kq + Kq~ decays, with the decay amplitudes shown in Eqs. fl30|) and 
( 1321) , we can find the numerical results in both scenarios for the CP- violating asymmetries 
in the pQCD approach follows, 



,dw n . v+v^ „ J -64.9^;°(o; 6 )^(/5)^; 9 7 (5f)^j(CKM)% (SI) 

5.9i^(c b )i^(/ 5 )il 1 9 2 ( J Bf)^i(CKM)% (S2) 



A***m ^k*+k*-\ ~ I - 49 - 9 -6:9(^)^5(/5)^: 8 7 (5f)^:?(CKM)% (si) 

' ' " 1 -98.9^M^(/ s )^( J Bf)^(CKM)% (S2) ■ (l>< 



a*ub ^k*+k*-) ~ / 12 - 9 ^:?(^)18:K/5)il:i(Bf)l8:I(CKM)% (si) 

CM 8 ° ° j ~ I -3.2l8;l( Wh8 )lg;i(/ s )i^(Bf )±£?(CKM)% (S2) ' m 



A"**m -> K*+K*-) ~ I 3.0ii;?(^)^(/ s )^(5f)^;J(CKM)% (SI) 
° " : ' " ~ 1 l-8^(^)^:i(/ S )t^(Sf)i°;{(CKM)% (S2) ' (,,) 



4 r« 7C+^ ~ / 99-liooK)loo(/5)l^W)l8J(CIKM)% (SI) 
AFsl S ° ° j ^ 1 99.9lSS( W6 .)jSS(/5)lSi(5f)lS5(CKM)% (S2) Ul ' 
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where the errors are induced by the uncertainties from shape parameter Ub(u)b s ) for B(B S ) 
meson, scalar decay constant f$ and Gegenbauer moments in the distribution ampli- 
tudes of scalar Kq, and CKM parameters (p,fj), respectively. For the direct CP asym- 
metries in the B d — > Kq + Kq~ and B s — > K£ + Kq~ decays for example, i.e., Eqs. ( |67|) 
and (1691) . one can find that their signs and magnitudes are rather different in both sce- 
narios within theoretical errors. In the former mode, A^p is (—69.4 ~ —57.6) in SI and 

(2.9 ~ 17.4) in S2; while in the latter one, Af P is (11.4 ~ 14.6) in SI and (-4.9 1.6) 

in S2. But, it is clear to find that the magnitudes of direct CP-violating asymmetries for 
these two decays in SI are much larger than those in S2 in the pQCD approach corre- 
spondingly, which will be more easily tested by the ongoing LHC and forthcoming SuperB 
experiments. Furthermore, once the predictions on Aq P including both size and sign in 
SI could be confirmed at the predicted level by the stringent experimental measurements 
in the future, which will be helpful to investigate the physical property of the scalar Kq 
and determine the better scenario describing its QCD dynamics in turn, and vice versa. 
Meanwhile, the stringent test of A% P {B d -> K^K^) can also provide indirect evidences 
for an important but controversial issue on the evaluation of annihilation contributions 
at leading power. 

Finally, it is worthy of mentioning that we here just study the perturbative short- 
distance contributions as the first step in present work. The large theoretical errors 
induced by the large uncertainties of the inputs in the nonperturbative distribution ampli- 
tudes, such as (pB q , 4>k*, etc, should be constrained by the precision measurements, which 
will be very helpful to explore the hadronic dynamics of Kq and the QCD dynamics in- 
volved in these considered decay channels. We do not consider the possible long-distance 
contributions, such as the rescattering effects, although they should be present, and they 
may be large and affect the theoretical predictions. It is beyond the scope of this work 
and expected to be studied in the future. 

V. SUMMARY 

In this work, we studied the two-body charmless hadronic B q — > K^K^q = u,d,s,c) 
decays by employing the pQCD approach based on the kr factorization theorem. Based 
on the assumption of two-quark structure of the light scalar we made the theoreti- 
cal predictions and phenomenological analysis on the physical observables: CP-averaged 
branching ratios and CP-violating asymmetries. 

From our numerical evaluations and phenomenological analysis, we found the following 
results: 

• In both scenarios of Kq, B q — > K^K* decays in the pQCD approach exhibit the 
large CP-averaged branching ratios (1CT 6 ~ 1CT 4 ), which are clearly measurable in 
the present B factories and ongoing LHC experiments. 

• In the considered six decay channels, only preliminary upper limit on Br(B^ — > 

Kq°K* ) m °de has been reported. The pQCD prediction basically agrees with this 
upper limit, and will be tested when better experimental measurements become 
available. 
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• For B u — > Kq + K decay, the large branching ratio plus large direct CP-violating 
asymmetry will offer a great opportunity to test the hadronic dynamics of Kq and 
the QCD dynamics involved in the considered mode. 

• The obvious nonzero A^p for B^s) — >■ Kq°K* decays will provide a good platform 
to explore the exotic new physics effects beyond the SM. 

• The pure annihilation process B c — > K$ K Q have a large branching ratio and will 
be measured soon at the LHC experiments, which can help us to understand the 
role of the annihilation contributions in B physics. 

• The pQCD predictions still have large theoretical uncertainties induced by the un- 
certainties of input parameters, such as the universal distribution amplitudes (f>B q 
and 4>k£, which should be well constrained by the precision data. 
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